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文章 编号 :1005-3085(2011)01-0101-08 
随机 微分 方程 解 的 二 次 型 估计 


马 洪 强 !， Rat 
(1- 复旦 大 学 上 海 视觉 艺术 学 院 ， 上 海 201620; 2- 东 华 大 学 应 用 数学 系 ， 上 海 201620) 
搞 要 : 由 于 随机 微分 方程 (SDE) 的 解析 解 求解 困难 ， 所 以 推导 SDE 解 的 不 等 式 估计 式 是 十 分 必要 的 . 
在 随机 系统 的 稳定 性 分 析 和 控制 设计 中 ， 李 亚 普 诺 夫 函数 常常 采用 二 次 型 函数 . 本 文 把 SDE 解 
的 传统 的 欧 几 里 德 范 数 形式 估计 式 推 广 到 SD 巨 解 的 二 次 型 估计 式 ， 包 括 解 的 矩 估计 和 几乎 必然 
估计 .我们 分 别 在 加 权 线 性 增长 条 件 和 加 权 单 边 增长 条 件 下 给 出 了 二 次 型 矩 估计 式 以 及 样本 李 亚 


普 诺 夫 指数 的 上 界 表 达 式 . 
关键 词 : 随机 微分 方程 ， 矩 估计 ; 几乎 必然 估计 ; 二 次 型 
分 类 号 : AMS(2000) 60H20 中 图 分 类 号 : 0211.63 文献 标识 码 : A 
1 引言 


为 了 研究 随机 微分 方程 (SDE) 的 稳定 性 ， 对 随机 微分 方程 的 解 进行 估计 其 有 重要 意义 . 随 
机 微分 方程 解 的 Euclid 范 数 形式 估计 已 被 广泛 研究 和 ， 但 是 对 有 关 随 机 微分 方程 解 的 其 它 形式 
的 估计 (比如 对 解 的 二 次 型 估计 ) 还 未 见报 道 . 

本 文 所 研究 的 问题 都 是 在 带 o 流 {Fi}:>o 的 完备 概率 空间 (N, F, P) 上 进行 的 ， 文 中 用 到 的 
标记 L?(Q, R") 表示 所 有 属于 OR” RAE EEP < 00 的 随机 变量 上 组 成 的 集合 ，.2p([a, b), R”) 
表示 所 有 满足 条 件 

/ | 的 Pdt < œ, as. 
Hn HEF, 适应 过 程 {j 归 }ja<i<s 组 成 的 集合 ，|z| 表 示 向 量 z 的 Eulid HR. 

令 丰 :BRxfto 如 一 Re 和 9:Rax[to7T] 一 Rxm Borel YR, B(t) An Eti BIZ). 
考虑 随机 微分 方程 


dx(t) = f (x(t), t)dt + g(a(t),t)dB(t), to <t<T, (1) 
z(to) = To. 


文献 [了 给 出 下 列 解 的 p 阶 和 矩 估计 和 几乎 必然 估计 . 
引 理 1 p>2, zo € L?(Q; R4)， 若 单 边 增长 条 件 成 立 ， 即 存在 一 个 常数 a > 0 使 得 对 所 有 
的 (zx,t) € R? x [to,T]， 有 


aT (1) fle, t) + >> *lol@,t)P < a( + e), (2) 
那么 
Ela(t)|? < 272" (1+ Elzo|?)e?o¢-), (3) 
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引 理 2 在 单 边 增长 条 件 (2) 下 ， 存 在 一 个 常数 a > 0， 随 机 微分 方程 解 的 样本 Lyapunov 指 
数 为 
lim sup > log |z(t)| < a. (4) 
too 


在 实际 应 用 中 ， 我 们 常常 要 考虑 SDE 解 的 函数 形式 ， 例如， 在 随机 稳定 性 分 析 和 控制 领 
域 ， 有 关 SDE 解 的 二 次 型 十 分 重要 中 ， 特 别 是 用 Lyapunov 直接 法 研究 稳定 性 时 ， 常 常 使 用 
二 次 型 &5， 而 线性 二 次 型 高 斯 (LQG) 问题 和 最 优 控制 问题 中 的 性 能 指标 函数 也 是 用 二 次 型 来 
表达 [le 站 ， 另 外 距离 判别 分 析 中 对 欧 氏 距离 加 权 后 构成 的 马 哈 拉 诺 比 斯 (Mahalanobis) 距离 也 是 
有 关 二 次 型 的 问题 由 等 等 . 

2 阶 矩 估计 和 几乎 必然 估计 在 研究 随机 z 阶 矩 指数 稳定 和 几乎 必然 指数 稳定 是 非常 重要 
的 . 在 随机 控制 领域 ， 用 Lyapunov 直 接 法 研究 稳定 性 上 F 蕊 ， 我 们 常常 用 二 次 型 V(z,t) = 

Z7(t)Qz(t) 来 代替 一 般 的 Lyapunov WATA, For @ 为 正定 和 矩阵， 那么 对 解 的 二 次 型 形式 
表达 式 的 性 质 分 析 就 显得 尤其 重要 . 本 文 在 文献 [1] 的 结果 基础 上 做 了 推广 ， 在 第 二 节 给 
出 SDB 解 的 二 次 型 p 阶 矩 估计 式 ， 在 第 三 节 证 明了 SDE 解 的 二 次 型 几乎 必然 估计 式 . 
2 p 阶 矩 估 计 


假设 z(t) 为 方程 (1) 的 解 ， 下 面 我 们 来 研究 关于 解 的 二 次 型 p 阶 矩 问题 . 
定理 1 p2 2, £o € L (Q; RÌ), 车 存在 一 常数 a， 则 对 所 有 的 (x,t) E€ Rd x [to, T]; A 


at’ Qf(a,t) + 2 二 trace g(x,t)" Qg(z,t) < a(ll + x7 Q7), (5) 
其 中 Q 为 对 称 非 负 定 矩阵 ，trace 表示 和 矩阵 的 秩 ， 那 么 对 任意 的 t> to A 
E (a7 (t)Qz(t))? < 2 (1 + E(x Qro)? Jeret- (6) 
证 明 WH, X (1+ 27 (t)Qax(t))? Alto 公式 得 
d(1+27(t)Qzx(t))? = p(1t+27(t)Qz(t)) > “al (t)Qf (a(t), t) 

+5 trace (9 (c(t), t) (plp — 2)(1 + 2? ()Qa(t)) "F(a" HQ)Pa" Q 
+p(1 +27 (t)Qu(t)) 五 Bee ))de 
+p(1+27(t)Qz(t)) Z “aT (t)Qg (x(t), t)dB(t) 

= p(1+27(1)Qa(t)) = ss t) 


E E = 2 |x? (HQg (e(t), #)? 
+5 trace (97 (a(t), DRIED) + aay ) 
+p(1+ 27 (t)Qzx(t))? “2T(t)Qg(z(t), t)dB(t). (7) 


由 于 存在 矩阵 U 使 得 
Q=U7U, 
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则 
(1+ 27 (t)Qz(t)) ‘|x? (t)Qg(a(t), t)? 
< (1+ aT (t)Qr(t)) Urt) PU galt), OF 
(1+ x7 (t)Qzx(t))~ (aT (t)Qzx(t)) trace g7 (2(t),t)Qg(x(t), t) 


< trace (97 (x(t), t)Qg(z(t), t)), 


lA 


A 


(1+ 27 (t)Qz(t))? 


< (1+ 0$Qro)¥ +p | (1+27(s)Qe(s))™ (TOREO, s) 





+P : trace g” (x(t), s)Qg(z(t), s) )ds 
; 2=2 
+p f (1+27(s)Qz(s)) ? x7(s)Qg(x(t), s)dB(s). 
to 
对 每 一 整数 n> 1， 定 义 停 时 
Tn =TAinf {t € [to, T] : |x| > n}. 
显然 7% TT, as. 由 条 件 (5) 和 It6 积分 性 质 ， 有 


E(1+27(t A Tp)Qx(t A t™m)) 3 


IA 


277 (1 + E3 Qzo)*) + Epa | "(1+ (a(s A Ta)" Qa(s A 7a))Ž)ds 
to 


= 273" (1+ E(xd Qao)?) + pa | E(1+ (27(s A m)Qz(s A Tn)) *)ds, 


由 Gronwall KÆRU, 4 


E(1+27 (tA t)Qzx(t A Tn)) 5 <27 (1+ E(t Qro)? )ertto). 
当 n 一 oo 时 ， 有 


B(e7 (t)Qu(t))® < E(1 +27 (t)Qa(t))® < 2° (1 + E(f Qro)? era". 


(8) 


(10) 


(11) 


推论 1 p> 2, zo € LPO; RI) KALRBR QAIEREIER, SDAA EKA 


(f7 (2, QF (@,t)) V (97 (z, t)Qg(z,t)) < KO +27Qz) 
成 立时 ， 则 有 不 等 式 (5) 成 立 ， 且 


a= VK + ECD. 
2 


(12) 
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证 明 当 条 件 (12) 成 立时 ， 式 (5) 中 的 
a=VK+ —— 
事实 上 ， 由 条 件 (12) 和 基本 不 等 式 2ab < a? +b, UHE > 0, FEMUR Q=UTU, K 


22T (t)Qf (a(t), t) < 2(VE waon (202 


A 


ven (EL 


< (ex? (t)U7Uc(t)) 4 ee 





< e(x7(t)Qz(t)) + eee (13) 
令 e = 二 VK， 得 
2 (Qf (a(t),t) < VE (1+27()Qz(t)), 
则 
aT (t)Qf(a(t),t) + P z T (x(t), t)Qg(x(t), t) 
< (VK + 0 (1+ 27 EQ). (14) 
3 “几乎 必然 估计 
考虑 随机 微分 方程 
dx(t) = f (z(t),t)dt + g(x(t),t)dB(t), te [to, 00), 
初 值 条 件 为 z(to) = zo € L?(0; R”). 假设 该 方程 在 [to, 00) 上 有 唯一 一 个 全 局 解 z(t). 
定义 1 依赖 于 二 次 型 的 样本 路 径 
lim sup > log (27 (t)Qz(t))*, 
称 为 加 权 样 本 Lyapunov 指数 ， 其 中 Q 为 非 负 定 矩 阵 . 
令 0<2Dp 入 2， 由 定理 1 可 知 ， 加 权 后 的 解 的 pz 阶 和 矩 满足 
E(27(t)Qz(t))? < P a + E(a% Qao)) Fert"), t> to. (15) 


这 就 意味 着 解 的 二 次 型 p 阶 矩 是 最 多 以 参数 为 pa 的 指数 增长 . 因此 ， 我 们 可 以 有 另 一 种 表达 
形式 
lim sup Z log (E(x? (t)Qz(t)) 3) < pa. 
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上 式 左 端 当 和 矩阵 Q@ ARMEEN, RIEZ A pE Lyapunov g% LRE H pE Lyapunov 
指数 不 会 超过 pa. 在 本 节 中 ， 我 们 将 做 解 的 几乎 必然 估计 . 首先 ， 我 们 估计 


lim sup > log (x7 (t)Qz(t)) 2 
t 一 Oo 
引 理 3 (指数 蒜 不 等 式 ) 令 g = (9 ,gn) € L?(R, R”), T, a, 8 是 任意 的 正 数 ， 那 么 


P{ sup | 了 9(s)dB(s) 一 ef la(s)Pas| > a} < eb, 


to<t<T o 


定理 2 在 加 权 单 边 增长 条 件 (5) 下 ， 随 机 微分 方程 解 的 加 权 样 本 Lyapunov 指数 将 不 会 大 
于 aw， 即 对 任意 的 上 之 如 ， 有 


lim sup = log (2T (t)Qzx(t))? <a, as. (16) 
证 明 由 It6 公式 和 加 权 单 边 增长 条 件 (5)， 有 
d(1+ (27 (t)Qz(t))) 
= (2x7 (t)Qf (a(t), t) + trace g7 (x(t), t)Qg(x(t), t))dt + 227 (t)Qg (z(t), t)dB(t). 


dlog (1 + (x7 (t)Qz(t))) 


= ( 227 (t)QF (x(t), t) + trace (97 (x(t), t)Qg(2(t), t)) 
1+ (27 ()Qz(t)) 





1 
—2 trace (Te), NR aO rga” Hemel, t)) )dt 
227 (t)Qg(x(t), t) 
Hrer) 2 ca 
两 边 积 分 ， 得 
log (1 + (x7 (t)Qz(t))) 


t 


= op (ite One) A Ta) CH Oee) r) 


. t | T 2 
+trace (g7 (z(7),7)Qg(z(7),7)))dr —2 ; 有 


”2Z7T(r)@o(z(r),r) 

+2 J ITOR PO 
t |zT(r)Qg(z(7), r)? 

< log (1+ (8 Quo) + 2a(t— to) -2 | TF et ar) 


t ah (r)Qa(2(r)-") anc, 
2 f IF AO EO =) 
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mo =2 f it? Oa or 


由 引 理 3， 对 任 一 整数 m > to 可 得 


|z7(r)Qg(x(r),7)/? 1 
PC me (ono -2 a > 26") <a 


由 Borel-Cantelli 定理 由 ， 对 几乎 所 有 的 w e Q， 存 在 一 个 随机 整数 no = now) > 如 十 1， 使 得 


t T r), 2 
2p, Mo -2 Tr ED] em nzn 


即 


i oe a to <t<n, n>no, a.s.. (19) 


因此 ， 对 几乎 所 有 的 we Q， 有 
3 log (1 + (x7 (t)Qz(t))) 


M(t) < 2logn+2 


IA 





1 
=z [log(l + (z6 Qzo)) + 2a(n — to) + 2logn], n>no, n-1<t<n, 


. 1 T 2 
lim sup z log (x7 (t)Qz(t)) 
t 一 OO 
1 
< limsup = log (1+ 27 (t)Qz(t)) 
toa at 


< limsup D [log(1 + 23 Qro) + 2a(n — to) + 2logn] =a, a.s.. (20) 
一 OO 


1 
2(n — 
考虑 加 权 线 性 增长 条 件 (12)， 存 在 一 常数 天 > 0 使 得 对 所 有 (x,t) © Rd x [to,00), A 

£7 ((, HQF (2, t)) V (97 (2, t)Qg(z,t)) < K(1+ 27Qz). 


易 知 ， 在 加 权 单 边 增长 条 件 中 a = VK +S. 因此 ， 我 们 容易 证 明 下 面 的 推论 . 

推论 2 在 加 权 线 性 增长 条 件 (12) 下 ， 方 程 (1) 的 解 的 加 权 样 本 Lyapunov 指数 不 超过 VK+ 
K, Bp 
2 


lim sup > log (a? (t)Qzx(t))? <VK+ > 
最 后 我 们 指出 ， 在 加 权 线 性 增长 条 件 (12) 下 ， 记 


2 (r)Qolz(t).r) 
my =2f ENOO 


为 连续 局 部 著 ， 那 么 平方 变 差 为 


t zT r T r 2 
(M(t), M(t)), =4 o ae < 4K (t — to). 
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因此 
lim sup (MO, MON <4K, as. 
zt 一 Oo 
则 
lim sup 地 一 0， .S.， 
t-00 t 
但 是 ， 加 权 单 边 增 长 条 件 并 不 能 保证 
lim sup MO =0, as 
too t 
成 立 . 
4 小 结 


本 文 是 在 随机 微分 方程 解 的 向 量 形式 的 基础 上 把 向 量 形式 推广 到 二 次 型 ， 把 对 解 


HY Euclid 范 数 形 式 的 估计 推广 到 二 次 型 ， 这 样 有 助 于 解决 控制 领域 中 的 稳定 性 问题 ， 使 得 
在 应 用 Lyapunov 直接 法 解决 稳定 性 问题 时 ， 用 二 次 型 来 代替 Lyapunov 函数 是 可 行 的 ， 这 样 就 
有 利于 找到 确定 的 Lyapunov 函数 ， 从 而 使 问题 具体 化 ， 也 方便 了 控制 工程 方面 人 员 的 研究 工 


作 ， 


这 也 是 本 文 的 意义 所 在 . 
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Quadratic Estimation to Solution of Stochastic Differential Equations 
MA Hong-qiang!, HU Liang-jian?t 


(1- Shanghai Institute of Visual Art, Fudan University, shanghai 201620; 
2- Department of Applied Mathematics, Donghua University, Shanghai 201620) 


Abstract: Since most stochastic differential equations (SDE) are not explicitly solvable, it is very 
important to find the estimation of the solution in the form of inequalities. In the research on stability 
analysis and control design of the stochastic systems, Lyapunov functions often take the quadratic 
forms. The aim of this paper is to extend the estimation from the classical Euclidean form to the 
quadratic form, including moment estimation and almost surely estimation of the SDE solution. As 
the results, the upper limits of moment estimation and sample Lyapunov index in quadratic function 
of solution are given under weighted linear growth condition and weighted one-side growth condition, 
respectively. 

Keywords: stochastic differential equations; moment estimation; almost surely estimation; quadratic 


form 
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